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Abstract Understanding the process of adaptation 
is a key mission in modern evolutionary biology. 
Animals living at high elevations face challenges in 
energy metabolism due to several environmental 
constraints (e.g, oxygen supply, food availability, 
and movement time). Animal behavioral processes 
are intimately related to energy metabolism, and 
therefore, behavioral modifications are expected 
to be an important mechanism for high-elevation 
adaptation. We tested this behavioral adaptation 
hypothesis using variations of motion visual 
displays in toad-headed agamid lizards of the genus 
Phrynocephalus. We predicted that complexity of 
visual motion displays would decrease with the 
increase of elevation, because motion visual displays 
are energetically costly. Displays of 12 Phrynocephalus 
species were collected with elevations ranging from 
sea level to 4600 m. We quantified display complexity 
using the number of display components, display 
duration, pathways of display components, as well 
as display speed for each species. Association between 
display complexity and elevation was analyzed using 
the phylogenetic generalized least squares (PGLS) 
model. We found that both the number of display 
components and the average value of tail coil speed 
were negatively correlated with elevation, suggesting 
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that toad-headed lizards living at high-elevation 
areas reduced their display complexity to cope with 
the environmental constraints. Our research provides 
direct evidence for high-elevation adaptation from a 
behavioral aspect and illustrates the potential impacts 
of environment heterogeneity on motion visual display 
diversification. 


Keywords high-elevation adaptation, lizard, motion visual 
display, Phrynocephalus, signal complexity 


1. Introduction 


Understanding the process of adaptation is one key mission 
in modern evolutionary biology and high-elevation 
adaptation provides a fertile research ground (Cheviron and 
Brumfield, 2012; Storz, 2021; Sun et al., 2018; Yang et al., 2012). 
Environmental constraints at high elevations, including low 
oxygen partial pressure, low ambient temperature, strong 
ultraviolet radiation, and great daily and seasonal fluctuations, 
pose serious challenges to most animals, particularly for 
ectotherms who are more sensitive to environmental changes 
(Hao et al., 2019). Nevertheless, many species inhabit high- 
elevation environments (Myers et al, 2000), and they often cope 
with these challenges by either suppressing total metabolism 
and oxygen demand, or increasing oxygen delivery efficiency 
and energy conversion (Li et al, 2018; Ramirez et al, 2007; Storz 
et al, 2010; Zhao et al, 2020). Animal behavioral processes have 
an intimate relationship with energy metabolism (Biro and 
Stamps, 2010; Kotiaho et al, 1998; Mowles, 2014; Ros et al., 2006; 
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Su et al., 2020), and animals at high-elevation areas likely adjust 
their behavior to adapt the environment by either constraining 
the activity intensity or shortening the activity time (Wu et 
al, 2018; Zhu et al, 2020). For example, Anolis lizards spent less 
time for activity and chose more open habitats with boulders 
for thermoregulation at high elevations (Muñoz and Losos, 
2017). Compared with physiological and genetic adaptation 
(Beall, 2007; Qu et al., 2021; Yang et al, 2014; Yang et al, 2016), 
behavioral adaptation at high-elevation environment has rarely 
been examined. 

As a ritual behavior, social displays transmit specific 
information during animal communication (Laidre and 
Johnstone, 2013). They play an important role in courtship, 
intra-sexual competition as well as territory defense, and 
their functions depend largely on complexity (Ord et al., 2001; 
Patricelli, 2016; Petrie et al., 1991; Wu et al, 2018). A diversity of 
organisms rely heavily on communication with a variety of 
different modalities (e.g, motion visual displays, olfactory, vocal; 
Endler, 1992; Laidre and Johnstone, 2013; Nelson and Jackson, 
2007). Physical movement-based motion visual displays are 
common in lizards’ communication (Bian et al., 2019; Fleishman 
and Pallus 2010; Rosenthal, 2007), and their complexity varies 
greatly among closely-related species due to different selection 
pressures (Nadhurou et al., 2016; Ord and Martins 2006). 
These motion visual display signals are energetically costly 
and complex displays often stimulate anaerobic metabolism 
(Bennett et al., 1981; Biro and Stamps, 2010; Clark, 2012; Ros 
et al., 2006; Wu et al., 2018; Zhu, 2020; Zhu et al., 2021). For 
example, compared with resting and moving, male wolf spiders 
(Hygrolycosa rubrofasciata) increase their metabolic rate 22 times 
during drumming behavior (Kotiaho et al., 1998). However, 
long anaerobic metabolism produces high concentrations of 
acid, which can affect blood and muscle pH, disrupt enzyme 
function and oxygen transport (Bennett, 1978). For high- 
elevation dwellers, we would expect them to reduce their 
motion visual display complexity to cope with the challenging 
environment and to avoid adverse impact from anaerobic 
metabolism. 

Quantifying display complexity can be difficult, and 
three methods are often used. The number of distinct displays 
or components in displaying repertoire is the most classical 
measurement, and has been widely used in assessment of 
display complexity in Anolis lizards (Decourcy and Jenssen, 1994; 
Greenberg and MacLean, 1978; Jenssen, 1978; Ord and Martins, 
2006;). For example, Martins et al. (2004) used the number of 
head-bobbing displays as an indicator of signal complexity. 
With advancement of the information theory, a system- 
based method has recently been used, in which the sequence 
variation among display components in succession was used 
as an indicator of display complexity (Fischer et al., 2017; Ord 
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and Martins, 2006; Shannon and Weaver, 1949). The higher 
the variation, the higher the signal complexity and the more 
information it may transmit (Peckre et al, 2019). Furthermore, 
variations within display components (e.g., speed, duration) 
have also been used as indicators of signal complexity (Freeberg 
et al, 2012; Hammerschmidt and Fischer, 2008; Vehrencamp, 
2000). 

The toad-headed agamid lizards (genus Phrynocephalus) 
provide an excellent model system to explore the display 
behavior adaptation to high-elevation environments. These 
lizards are widely distributed in a large elevational range from 
sea level to 4 600 m and display a variety of motion visual 
signals. For example, both males and females of P. vlangalii, a 
high-elevation species, use tail coil and tail lash in territory 
defense and courtship, and the speed of tail display likely 
encodes levels of individual threat (Qi et al., 2011; Qiu et al., 2018). 
Several lowland Phr ynocephalus species show more complex 
displays besides tail coil and tail lash, such as limbs flapping 
and body turning (Lin, 2020) In addition, the athletic ability of 
P. vlangalii depends largely on elevation, and individuals from 
high altitudes have relatively weak locomotor performance 
compared with low-elevation individuals (Wu et al., 2018a). 
Furthermore, our recent studies have shown that speed of 
tail display in P. vlangalii is mainly regulated by anaerobic 
metabolism, and fast display speed is associated with high blood 
lactate concentration (Zhu et al., 2020; Zhu et al., 2021). 

In this study, we explore the relationship between high- 
elevation adaptation and complexity of motion visual display 
in the genus Phrynocephalus. We collected display data from 
12 species and quantified the display complexity using the 
motion visual signal digitization method (Hedrick, 2008; Peters 
et al., 2016). The association between display complexity and 
elevation was analyzed using the phylogenetic generalized least 
squares (PGLS) model. We predict that complexity of motion 
visual displays would decrease with elevation, because motion 
visual displays are energetically costly and total metabolism is 
constrained at high elevations. 


2. Materials and methods 


2.1. Species and study sites We collected motion visual display 
data from 12 Phr ynocephalus species in the Gobi Desert of 
northwestern China and the Qinghai-Tibetan Plateau (Figure 
1) Lizards in different social context produce different signal, 
and all displays in this study were under a male-male social 
context, i.e, male residents competing with male intruders. 
Habitats for each species are similar and mostly arid, varied 
from rock substrate to sand dunes, with sparse and low 
vegetation. Elevation of study sites varied from 604 m (P. 
alpherakii) to 4 550 m (P. er ythrurus). For two widespread species, 
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P. axillaris and P. versicolor, we examined two populations and 
their elevations were averaged between the populations. The 
study site and sampling information is presented in Table 1. 


22. Display signal data collection We first located and sexed 
target lizards using binoculars. Once a target was selected, one 
investigator (QH) set up a video camera (Sony HDR PJ670) 
at approximately three meters away from the target. After a 
five minute acclimation period, a second investigator (XQ) 
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introduced an intruder lizard toward the target from four 
meters away. The intruder was tied with a 30-centimeter-long 
dental floss around the waist and was tethered with a 4-meter 
fishing rod, which allowed the lizard to move freely. At the 
same time, the first investigator started the camera and filmed 
the motion visual displays of the target lizard. We would end 
the trial if there was a potential conflict escalation or target 


did not show displays within two minutes. Immediately after 


100°E 


Figure 1 Map of study sites and the image of 12 Phrynocephalus species examined in this study. 1. P. mystaceus, elevation 604 m. 2. P. 
alpherakii, elevation 604 m. 3. P. helioscopus, elevation 637 m. 4. P. forsythii, elevation 894 m. 5a. P. axillaris, elevation 894 m. 5b. P. axillaris, 
elevation 1 150 m. 6a. P. versicolor, elevation 870 m. 6b. P. versicolor, elevation 1 872 m. 7. P. przewalskii, elevation 1 393 m. 8. P. putjatai, 
elevation 2 303 m. 9. P. guinanensis, elevation 3 019 m. 10. P. vlangalii, elevation 3 457 m. 11. P. erythrurus, elevation 4 550 m. 12. P. 


theobaldi, elevation 3 670 m. 


Table 1 Study site and sample size information of the 12 Phrynocephalus species examined in this study. 


Id Species Location Coordinates Elevation (m) Sample size 

1 Phrynocephalus mystaceus Xinjiang Uygur Autonomous Region, China 80.8131°E, 44.0095°N 604 24 
Phrynocephalus alpherakii Xinjiang Uygur Autonomous Region, China 80.8131°E, 44.0095°N 604 15 
Phrynocephalus helioscopus Xinjiang Uygur Autonomous Region, China 83.5930°E, 46.1912°N 637 24 
Phrynocephalus forsythii Xinjiang Uygur Autonomous Region, China 86.2239°E, 41.4975°N 894 17 

5a Phrynocephalus axillaris Xinjiang Uygur Autonomous Region, China 88.0275" E, 38.9713°N 894 9 

5b Phrynocephalus axillaris Gansu Province, China 94.6567" E, 40.1004°N 1150 8 
6a Phrynocephalus versicolor Xinjiang Uygur Autonomous Region, China 93.6219°E, 43.0453°N 870 iu 
6b Phrynocephalus versicolor Gansu Province, China 95.4506°E, 41.1074°N 1 872 12 
Phrynocephalus przewalskii Gansu Province, China 102.9817E, 38.4950°N 186595; 15 
Phrynocephalus putjatai Qinghai Province, China 101.2419" E, 36.0797N 2 303 10 
Phrynocephalus guinanensis Qinghai Province, China 100.6939°E, 36.0797°N 3019 16 
10 Phrynocephalus vlangalii Zoige County, Sichuan Province, China 102.9213'E, 32.9333°N 3 457 39 
11 Phrynocephalus erythrurus Qinghai Province, China 92.4853°E, 34.2336 N 4550 34 
12 Phrynocephalus theobaldi Xizang (Tibet) Autonomous Region, China 90.9888 E, 29.5804" N 3670 27 
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filming, we captured the target lizards and measured their body 
temperature. To alleviate the effect of individual physiological 
condition on display complexity, we confined our display 
trials only during 11 am to 4 pm in sunny days when lizards 
are most active. To avoid impacts from potential previous 
social interactions and individual body size, a size-matched 
intruder collected from a different population was used for 
each target. After each trial, a ping-pong ball was placed at the 
exact location of the target and was filmed to serve as a scale in 
subsequent display digitization. 


2.3. Display complexity measurement We measured 
display complexity for each species using five parameters: the 
number of display components, display duration, number of 
pathways among different components, as well as tail coil 
and tail lash speed (Ord and Martins, 2006; Shannon, 1948). 
The component was defined as specific and repeatable display 
posture (Ramos and Peters, 2017). Some components were 
shared by several species and others were species-specific, and 
therefore, the number of components represented an important 
aspect of display complexity (Freeberg et al., 2012; Miller and 
Osmanski, 2009; Vehrencamp, 2000). The display duration 
was defined as the average time sustained by a species-typical 
display (Hammerschmidt and Fischer, 2008; McComb and 
Semple, 2005). For each species, the display sequence was mostly 
fixed, but individuals sometimes omitted some components 
under specific social context (Fischer et al., 2017). To establish 
the species-typical display sequence, we analyzed the motion 
visual display bout by bout and determined the starting 
and terminating component respectively for each species. 
We then calculated the transition probability (probability 
of one component followed by the other) between different 
components, and the display sequence was described as the 
sequence with the highest transition probability among 
components. The number of pathways between different 
components was determined using the number of different 
transition probability. The average speed of tail coil and tail lash 
were quantified following the methods outlined by Hedrick 
(2008) and Peters et al. (2016). Briefly, we extracted video footages 
of tail displays using the program iskysoft (iSkysoft Technology 
Corp.) and tracked the movement of tail tips in Matlab 2016b 
(MathWorks Inc, Natick, MA, USA). The position of the tail 
tip was located in successive frames to generate x-y coordinate 
data over time for each display. The x-y coordinates were then 
converted to millimeter using the ping-pang ball in the image as 
a scale (Peters et al., 2016; Wu et al., 2018). The Euclidean distance 
between successive digitized position points provided a vector 
of speed measurements for the whole sequence. We averaged 
these across the display sequence to determine the average 
speed. To reflect the potential difference between components, 
we calculated the speed for tail coil and tail lash respectively. 
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To account for the effect of orientation of lizard relative to the 
camera, we categorized each display as either facing towards/ 
away from the camera or at right angles to the camera, 
and calculated the display speed from different orientation 
respectively (Bian et al., 2016). 


2.4. Data analysis To examine the association between display 
complexity and elevation, we tested the relationships of the 
number of display components, display duration, number of 
pathways among different components, tail coil speed, and tail 
lash speed against elevation. To account for the phylogenetic 
non-independent effect, we used the phylogenetic generalized 
least-squares (PGLS) model in package “caper” (Orme et al., 2013). 
Elevation was treated as the predictor variable, while number 
of components, display duration, number of pathways among 
different components, tail coil speed, and tail lash speed were 
treated as the response variable, respectively. A phylogenetic 
tree of the 12 study species (Figure 2) was constructed and 
modified from Solovyeva et al. (2018). We only used the display 
speed from facing towards/away from the camera to avoid the 
effect of orientation of lizards. 

To test the phylogenetic effect on display complexity, we 
also estimated Pagel's 4 for the number of display components, 
display duration, number of component connections, tail coil 
speed, and tail lash speed. Pagel's 2 was estimated using the 
“phytools” package (Lynch, 1991; Pagel, 1999; Revell, 2012), which 
ranges from 0 to 1, with “O” meaning weak phylogenetic effect 
and “I” meaning strong phylogenetic effect (Freckleton et al., 
2015). Previous studies suggested that Pagel’s 4 is an effective 
measurement and typically performs better than other 
commonly used metrics for discriminating between more 
complex models of trait evolution (Miinkemiiller et al., 2012). All 
statistical analysis were conducted in R 362 (R Development 
Core Team, 2019). 


2.5. Ethical approval All applicable international, national, 
and institutional laws and guidelines for the care and use 
of animals were strictly followed. All activities were under 
permission from local conservation authorities and animal 
handling followed the approved protocols (protocol number 
2017005, Chengdu Institute of Biology). 


3. Results 


3.1. Signal repertoire and species-typical sequence We 
analyzed a total of 261 display bouts from 131 individuals of 
12 species. Among all those displays, 16 different components 
were defined, including tail lashing (TL), tail coiling (TC), “8” 
tail lashing (8), standing up (SU), push-up (PU), leaning (LE), 
tail lashing with four limbs flapping (T+F), and several others. 
Table 2 presents a complete list of all components and their 
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Figure 2 A phylogenetic tree of the 12 Phrynocephalus species examined in this study (modified from Solovyeva et al., 2018), along with 
schematic diagrams of the component network graph and display components for each species. 
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Table 2 Sixteen visual display components of the 12 Phrynocephalus species examined in this study. 


Display component (Abbreviation) Definition Component diagram 


The upper limb elbow bends before the first half of the body falls and then 


Pushing up (PU) holds up 


Four limbs ground flapping (FGF) Four limbs flapping the ground alternatively while the body circling Sa 


“8” tail lashing with high coiling (S8U) Tail is raised above the height of the head and lashing in a “8” track 


Tail coiling (TC) Tail raises above the height of the head and curls 


Tail coils at a low place (LTC) Tail coils close to the ground 


Tail lashes at a low place (LTL) Tail lifts a little and close to the ground with lashing slightly 


Tail raises above the height of the head in the process of accompanying 


Tail raises with lashing (LU) låsbing at thesignal beginning 


Tail lashes + Four limbs ground Tail lashes above the height of the head with four limbs flapping the 
flapping (T+F) ground 
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definitions. Among them, the first three types (TL, TC, and 
8) were basic components, and they might give rise to other 
components. For example, tail lashing (TL) might appear as tail 
lashing at a low place (LTL), tail raising with waggling (HTW), 
tail raising with lashing (LU), or tail falling with lashing (LTD). 
Tail coiling (TC) might appear as tail coiling intermittently 
(IC) or tail coiling at a low place (LTC), while “8” tail lashing (8) 
might appear as “8” tail lashing with high coiling (S8U) or “8” 
tail lashing with low coiling (L8) (Table SD. 

Most species possess their own display repertoire and have 
a species-typical sequence. We provide a detailed description in 
Table S2. The index of display complexity, including number 
of components, display duration, number of pathways, average 
values for tail coil and tail lash speed are presented in Table 3. 


3.2. Correlation between display complexity and elevation 
We found significant and strong phylogenetic signal in 
display duration (A = 0.9918, P = 0.0090) and number of display 
components (A = 0.6000, P = 0.0266). The average speed of 
tail lash had a large å but a marginally significant P (A = 
0.9711, P = 0.0906), which we interpreted as having less strong 
phylogenetic signal. The average speed of tail coil had a large 4 
but insignificant P value (A = 0.8356, P = 0.9393), which could be 
due to the small number of species. The number of pathways 
among display components had a small 4 and an insignificant P 
value (A = 0.1872, P = 0.4771), suggesting minimum phylogenetic 
signal for this trait (Figure SI). 

For the association between display complexity and 
elevation, the number of display components (f = -0.0007, P = 
0.0063, A = 0) and average tail coil speed (8 = —0.0026, P = 0.0459, 
À =1) were negatively correlated with altitudes (Figure 3, Table 
4). We found no association of the display duration (f = —0.0016, 
P = 0.1985, A = 1), the number of pathways (f = —0.0008, P = 
0.2596, A = 0), the tail lash speed (8 = -0.0071, P = 0.1849, 2 = 1) 
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with altitudes (Table 4). 
4. Discussion 


Our results clearly demonstrated that Phrynocephalus species 
at high-elevation areas have modified their motion visual 
display signals. First, species at high elevations appear to 
have fewer display components in general (Figures 2-3). The 
negative correlation between elevation and the number of 
display components is clear and significant. Of all the high- 
elevation species, P. guinanensis has the fewest components, but 
it is not the highest species. This species has an obvious colorful 
belly and trunk side for both male and female individuals, 
which most other species do not have. We postulate that the 
presence of colorful belly may compensate for the reduction 
of display components. P. forsythii presents another interesting 
case. It is a part of the high-elevation clade but currently has a 
primary lowland distribution. It has a large number of display 
components similar to other lowland species. Its phylogenetic 
position suggests that it originated from a high-elevation 
ancestor and turned to lowland secondarily (Figure 2). Its large 
number of display components further supports the association 
between reduction in display component and high-elevation 
environment. Second, species at high elevations appear to 
have slower tail coil speed. Two other parameters, the display 
duration and tail lash speed, are also lower at high elevations, 
although the differences are not statistically significant. We 
suggest these modifications represent adaptations to high- 
elevation environments of these lizards. Other behavioral 
modifications were also reported in lizards. For example, the 
horned lizards (Phrynosoma hernandesi) adjust their basking 
duration to cope with low temperature challenges when 
transplanted to high-elevation areas (Refsnider et al., 2018), 
and Anolis lizards chose more open habitats with boulders for 


Table 3 Number of display components, average display duration (s), number of pathways among different components, and average values 
(cm/s) for tail coil speed and tail lash speed for each Phrynocephalus species. 


Species Number of Signal duration Number of Average speed of Average speed of 
components (+ SE) pathways tail coil (+ SE) tail lash (+ SE) 
P. mystaceus 7 49.48 +14.67 s 13 DOMI 38.66 13.23 
P. alpherakii 9 29.47 + 4.52 s 18 11.28 + 3.03 63.03+ 6.20 
P. helioscepus 8 TESI a8 ESS 15 16.66 + 0.70 16.44 € 1.43 
P. forsythii 7 20.96 + 2.41 s 13 21.08 + 1.73 28.05 + 1.68 
P. axillaris 8 18.99 + 3.[21s 13 Sakai) ae BL) 143.40 + 9.40 
P. versicolor 7 14.33 + 1.64s 18 10.76 + 1.34 63.37 + 4.93 
P. przewalskii 8 20.82 + 3.07 s 15 Wi ae 115227 70.42 + 10.28 
P. putjatia 6 12.66 + 1.85 s 12 15.59 + 1.17 27.19 + 3.44 
P. guinanensis 4 8.89 + 0.67 s 5) DIS 24.42 € 3.65 
P. vlangalii 6 15.91 + 3.14 s 16 27.77 t 2.76 12.92 + 1.79 
P. erythrurus 6 17.90 + 3.10s 3) 12572077 45.72 + 2.44 
P. theobaldi 6 16.46 + 1.85 s 13 12.39 £1.21 1.26 + NA 
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Figure 3 A: Correlation between number of display components and elevation (R° = 0.4961, P = 0.0063). B: Correlation between aver- 


age tail coil speed and elevation (R^ = 0.2760, P = 0.0459). 


Table 4 The association between display complexity and elevation across 12 Phrynocephalus species. The display complexity was measured 
using number of display components, display duration, and number of pathways. Significant predictors are marked in bold. 


2 


A b t R P 
Number of signal components 0 —0.0007 -3.4395 0.4961 0.0063 
Display duration 1 -0.0016 -1.3773 0.0754 0.1985 
Number of pathways 0 —0.0008 -1.95 0.0375 0.2596 
Average speed of tail coil 1 -0.0026 -2.2791 0.276 0.0459 
Average speed of tail lash 1 -0.0071 -1.4241 0.0855 0.1849 


thermoregulation at high elevations (Muñoz and Losos, 2017). 
The reduced visual display complexity at high elevations 
in Phrynocephalus lizards is likely caused by metabolic 
capacity constraints. Many ectothermic species living in harsh 
high-elevation areas have reduced metabolic capacity (e.g., 
reducing oxygen consumption and metabolic rate) due to 
environmental constraints (Tan et al. 2021). For example, the 
lactate dehydrogenase (LDH) activity, an indicator of anaerobic 
metabolism, is lower in P. erythrurus (elevation 5 300 m) than 
in P. przewalskii (elevation from 1 000 to 1 500 m) (Tang et al., 
2013). Furthermore, recent research showed that the intensity of 
tail displays in P. vlangalii (altitude from 2 000 to 4 500 m) was 
primarily regulated by anaerobic metabolism (Zhu et al., 2021). 
Therefore, reduction in their number display components and 
speed of display by the high-elevation Phrynocephalus species is 
most likely response to their low anaerobic metabolic capacity. 
Alternatively, energy constraints on the nervous system 
at high-elevation environments may also contribute to 
the simplification of visual displays. The central nervous 
system regulates the movement of skeletal muscles in most 
vertebrates (Akinrodoye and Lui, 2021; Kingsbury et al., 2019), 
and as accurate and complex communication signals, the tail 
displays of Phrynocephalus must be closely controlled by the 
brain. Nevertheless, the function of brain depends largely on 


continuous supply of blood sugar and oxygen (Ereciriska and 
Silver, 2001; Olesen, 1986). High-altitude hypoxia undoubtedly 
imposes various constraints on brain functions and oxygen 
conditions during development can affect lizards’ cognition 
as well (Sun et al., 2014), which may lead to reduced display 
complexity. This has been partially evidenced by recent works 
on the Asiatic Toad (Bufo gargarizans) that individuals at high- 
elevation areas have reduced brain sizes and limited movement 
capacity (Mai et al., 2017; Yao et al., 2020). Future research on 
relationships between brain size and display complexity, as well 
as detection of the critical functional brain regions associated 
with display manipulation is needed. 

One confounding issue with the observed motion visual 
display simplification in high-elevation Phr ynocephalus species is 
social complexity. According to the social complexity hypothesis 
(SCH), animals living in complex groups likely encounter more 
social interactions, and may evolve complex communication 
signals to cope with social conflict and maintain social 
coalitions (Freeberg, 2006; Freeberg et al., 2012). Our field 
observations show that Phrynocephalus species at high-elevation 
areas are more territorial (eg. stay near burrow and defend 
territory intensively) compared with species in lowland; based 
on SCH, we would expect that high-elevation species possess 


more complex signals. This is the opposite of our results in this 
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study. We assumed that this is likely a compensation way for 
the effect of display simplification, species from high-elevation 
area likely increase their territoriality and avoid unnecessary 
social conflict. Potential associations among social complexity, 
territorial defense, and display complexity remain to be 
explored and tested. 


5. Conclusion 


Lizards of the genus Phrynocephalus use complex and species- 
specific motion visual displays during social communication. 
Species living in high-elevation areas reduce their display 
complexity, particularly the number of display components and 
tail coil speed. This is likely associated with limited metabolic 
capacity and brain function. The genus Phrynocephalus provides 
an excellent system for studying high-elevation adaptation, and 
more research should be conducted on this group. 
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Figure S1 Phylogenetic signal in A: display duration; B: number of display components; and C: average tail lash speed. 


Table S1 The frequency of different components which is used as starting and terminating components and the transition probabilities 
between components for each Phrynocephalus species in China. Referring to Table 2 for the meaning of component abbreviation. 
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Table S2 The detailed description of species-typical sequence for each Phrynocephalus species. The components in parentheses of the typical 
sequences are the repeatable sections, and "/" represents "or". 


Species Typicalsequence Description 
P. mystaceus exhibited a typical display sequence of “LE-LU-8-(TC-LTD)-TC/TL/LTD”. Their display 
typically started with tilting the body left or right and showing the dorsum to the opponent (LE 79.17%). 
Then their tails raised above the height of the head in the process of accompanying lashing at the signal 
beginning (LU 12.30%), followed by tail raising above the height of the head in the process of accompanying 
LE-LU-8-(TC- lashing at the signal beginning (“8” 9.84%). Next, the tails raised above the height of the head and curling (TC 
LTD)-TC/TL/LTD 11.48%), followed by tail putting down from coil accompanying lashing (LTD 27.05%) or tail raising high 
above the height of the head and lashing sharply toward both body sides (TL 8.20%). They might display "TC” 
followed by “LTD” repeatedly. Their terminal action was tail putting down from coil accompanying lashing 
(LTD 37.50%) or tail raising above the height of the head and curling (TC 29.17%), or tail raising high above 
the height of the head and lashing sharply toward both body sides (TL 25.00%). . 


P. alpherakii exhibited a typical display sequence of “LU-(S8U-T+F/TL)-TC/T+F”. Their display typically 
started with tail raising above the height of the head in the process of accompanying lashing at the signal 
beginning (LU 73.33%) or tail lashing above the height of the head with four limbs flapping the ground (T+F 
26.67%). Then tail is raised above the height of the head and lashing in a “8” track (S8U 8.33%), followed 
by tail lashing above the height of the head with four limbs flapping the ground (T+F 21.21%), or just tail 
raising high above the height of the head and lashing sharply toward both body sides (TL 17.42%). They 
could display “S8U” followed by “T+F” or “TL” repeatedly. Finally, their terminal action was tail raising 
above the height of the head and curling (TC 53.33%) or tail lashing above the height of the head with four 
limbs flapping the ground (T+F 20.00%). 


P. helioscopus exhibited a typical display sequence of “SU-LU-IC-HTW-(TC-LTD)-LTD”. Their display 
typically started with walking forward and holding up body before signaling (SU 75.00%). Then the tail 
raising above the height of the head in the process of accompanying lashing at the signal beginning (LU 

SU-LU-IC-HTW- 3.92%), followed by tail raising and coiling intermittently (IC 3.92%), and then tail raising high above the 
(TC-LTD)-LTD height of the head and the tip of the tail waggling on both sides (HTW 7.84%), followed by tail raising above 
the height of the head and coiling (TC 9.80%), and then tail putting down from coil accompanying lashing 
(LTD 35.78%). They could display "TC” followed by “LTD” repeatedly. Their terminal action was tail putting 
down from coil accompanying lashing (LTD 70.83%). 


P. forsythii exhibited a typical display sequence of “PU-LU-IC-(HTW-TC)-TC”. Their display typically 
started with The upper limb elbow bends before the first half of the body falls and then holds up (PU 
58.82%). Then the tail raising above the height of the head in the process of accompanying lashing at the 
PU-LU-IC-(HTW- signal beginning (LU 8.82%), followed by the tail raising and coiling intermittently (IC 10.78%). And then the 
TC)-TC tail raising high above the height of the head and the tip of the tail waggling on both sides (HTW 21.57%), 
followed by the tail raising above the height of the head and curling (TC 23.53%). They could display “HTW” 
followed by "TC" repeatedly. Their terminal action was tail raising above the height of the head and curling 
(TC 82.35%). 


P. axillaris exhibited a typical display sequence of “LU-(S8U-T+F/TL)-TC/TL”. Their display typically 

started with tail raising above the height of the head in the process of accompanying lashing at the signal 

beginning (LU 35.29%) or started by tail is raised above the height of the head and lashing in a “8” track 
(S8U 29.41%), they also could walk slowly forward and hold up body before signaling (SU 29.41%). Then tail 

LU-(S8U-T+F/TL)- raising above the height of the head and lashing in a “8” track (S8U 9.09%), followed by tail lashing above 
TOANE the height of the head with four limbs flapping the ground (T+F 24.7%), or tail raising high above the height 
of the head and lashing sharply toward both body sides (TL 16.67%). They could display “S8U” followed by 

“T+F” or “TL” repeatedly. Finally, their terminal action was tail raising above the height of the head and 
curling (TC 76.47%) or tail raising high above the height of the head and lashing sharply toward both body 

sides (TL 17.65%). 


P. versicolor exhibited a typical display sequence of “LU-(S8U-T+F/TL)-TC”. Their display typically started 
with tail raising above the height of the head in the process of accompanying lashing at the signal beginning 
(LU 52.17%). Then their tails raised above the height of the head and lashing in a “8” track (S8U 7.43%), 
followed by tail lashing above the height of the head with four limbs flapping the ground (T+F 23.65%), or 
just tail raising high above the height of the head and lashing sharply toward both body sides (TL 15.54%). 
They might display “S8U” followed by “T+F” or “TL” repeatedly. Their terminal action was tail raising above 
the height of the head and curling (TC 73.91%). 


P. przewalskii exhibited a typical display sequence of “LU-(S8U-T+F/TL)-TC”. Their display typically started 
with tail raising above the height of the head in the process of accompanying lashing at the signal beginning 
(LU 66.67%) or tail raising above the height of the head and curling (TC 20.00%). Then tail raising above the 
height of the head and lashing in a “8” track (S8U 10.3456), followed by tail lashing with four limbs flapping 
the ground (T+F 15.52%), or just tail raising high above the height of the head and lashing sharply toward 
both body sides (TL 24.14%). They could display “S8U” followed by "T+F” or "TL" repeatedly. Their terminal 
action was tail raising above the height of the head and curling (TC 86.6796). 


P. mystaceus 


LU-(S8U-T+F/TL) 


P. alpherakii -TC/ T4F 


P. helioscepus 


P. forsythii 


P. axillaris 


LU-(S8U-T+E/TL) - 
TC 


P. versicolor 


LU-(S8U-T+F/TL)- 
TC 


P. przewalskii 


Continued Table S2 


Species Typicalsequence Description 


P. putjatia exhibited a typical display sequence of “LU-TC-HTW/TL-TC”. Their display typically started 
with tail raising above the height of the head in the process of accompanying lashing at the signal 
beginning (LU 50.00%) or walking forward and holding up body before signaling (SU 40.00%). Then the tail 
raising above the height of the head and curling (TC 20.6996), followed by the tail raising high above the 
height of the head and the tip of the tail waggling on both sides (HTW 13.7996) or tail raising high above 
the height of the head and lashing sharply toward both body sides (TL 10.34%). Their terminal action was 
tail raising above the height of the head and curling (TC 100.0056). 


LU-TC-HTW/TL- 


P. putjatia TC 


P. guinanensis exhibited a typical display sequence of “SU/TC-HTW-TC”. Their display typically started 
with walking forward and holding up body before signaling (SU 50.00%) or tail raising above the height of 
P. guinanensis | SU/TC-HTW-TC the head and curling (TC 50%). Then the tail raising high above the height of the head and the tip of the tail 
waggling on both sides (HTW 37.5096). Their terminal action was tail raising above the height of the head 
and curling (TC 93.7596). 


P. vlangalii exhibited a typical display sequence of "SU-LU-(TC-TL)-TC'". Their display typically started 
with walking forward and holding up body before signaling (SU 64.1096) or tail raising above the height 
of the head and curling (TC 23.08%). Then tail raising above the height of the head in the process of 
P. vlangalii ^ SU-LU-(TC-TL)-TC accompanying lashing at the signal beginning (LU 5.36%), followed by the tail raising above the height of 
the head and curling (TC 4.46%), and then tail raising high above the height of the head and lashing sharply 
toward both body sides (TL 25.8996). They could display "TC" followed by "TL" repeatedly. Their terminal 
action was tail raising above the height of the head and curling (TC 100.056). 


P. erythrurus exhibited a typical display sequence of "SU-LU-(TC-TL)-TC'". Their display typically started 
with walking forward and holding up body before signaling (SU 67.6596). Then tail raising above the height 
of the head in the process of accompanying lashing at the signal beginning (LU 12.8796), followed by the tail 
raising above the height of the head and curling (TC 15.84%), and then tail raising high above the height of 
the head and lashing sharply toward both body sides (TL 29.7096). They could display "TC" followed by "TL" 

repeatedly. Their terminal action was tail raising above the height of the head and curling (TC 85.2996). 


P. erythrurus SU-LU-(TC-TL)-TC 


P. theobaldi exhibited a typical display sequence of “TC/SU/LU-TC-HTW-TC”. Their display typically 
started with walking forward and holding up body before signaling (SU 37.04%) or tail raising above the 
TC/SU/LU-TC- height of the head in the process of accompanying lashing at the signal beginning (LU 22.22%), even tail 
P. theobaldi HTW-TC raising above the height of the head and curling all the time (TC 40.74%). Then tail raising high above the 
height of the head and the tip of the tail waggling on both sides (HTW 20.7556), followed by the tail raising 
above the height of the head and curling (TC 18.87%). Their terminal action was tail raising above the 
height of the head and curling (TC 100.0096). 


